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Abstract 

The influence of the chiral mean field on the K + transverse flow in heavy ion 
collisions at SIS energy is investigated within covariant kaon dynamics. For 
the kaon mesons inside the nuclear medium a quasi-particle picture including 
scalar and vector fields is adopted and compared to the standard treatment 
with a static potential. It is confirmed that a Lorentz force from spatial 
component of the vector field provides an important contribution to the in¬ 
medium kaon dynamics and strongly counterbalances the influence of the 
vector potential on the K + in-plane flow. The FOPI data can be reasonably 
described using in-medium kaon potentials based on effective chiral models. 
The information on the in-medium K + potential extracted from kaon flow is 

consistent with the knowledge from other sources. 
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I. INTRODUCTION 


Properties of kaons in dense hadronic matter are important for a better understanding of 
both, a possible restoration of chiral symmetry in dense hadronic matter and the properties 
of nuclear matter at high densities. It is known from chiral models that the kaon mean field 
is related to chiral symmetry breaking [1], The in-medium effects give rise to an attractive 
scalar potential inside the nuclear medium which is in mean field approximation proportional 
to the kaon-nucleon Sigma term Ekn- A second part of the mean field originates from the 
interaction with vector mesons [1-5]. The vector potential is repulsive for kaons K + and, 
due to G-parity conservation, attractive for antikaons K~. A strong attractive potential for 
antikaons may also favor K~ condensation at high nuclear densities and thus modifies the 
properties of neutron stars [6]. 

There have been extensive experimental efforts to search for the kaon in-medium prop¬ 
erties in heavy ion collisions, in particular at intermediate energies [7-11]. Corresponding 
transport calculations revealed significant evidence for in-medium modifications of the kaon 
properties during the course of such reactions [12-21], The picture was recently comple¬ 
mented by the measurements of K + production in proton-nucleus reactions at COSY-ANKE 
[22]. While A + A reactions test the kaon dynamics also at supra-normal nuclear densities, 
p + A reactions can not exceed saturation density. However, in the latter case the reaction 
dynamics is theoretically much easier to handle (clue to much less secondary scatterings) 
and therefore the interpretation of the data is more straightforward. From the p + A data 
the existence of a repulsive in-medium K + potential of about 20 MeV at saturation density 
(U K (po) — 20 MeV) was derived [22], 

In A + A reactions the transverse flow of K + mesons is one of specially attractive observ¬ 
ables. It was proposed by Ko and Li [12] that the kaon flow pattern at the final stage can 
be used as a sensitive probe for the kaon potential in a nuclear medium. Other theoretical 
studies predicted similar features for the kaon flow [14,18]. As pointed out by Fuchs at al. 
[19], however, the kaon dynamics used in the above investigations is non-covariant, i.e. it 


2 



is based on a static potential like force. The Lorentz-force like contribution from a typi¬ 
cal relativistic scalar-vector type structure is missing. However, this contribution strongly 
counterbalances the influence of the vector potential on the K + in-plane flow which makes 
it more difficult to draw definite conclusions from transverse flow pattern. 

In the meantime the FOPI collaboration [23] published new data of K + side flow with 
improved precision. These open the possibility to study this problem in more detail. In 
this work we apply the covariant kaon dynamics to describe the kaon in-medium properties, 
discuss its influence on the K + in-plane flow in heavy ion collisions, and compare to the 
standard treatment with a static potential. Our studies show that a relatively strong repul¬ 
sive K + potential as proposed by Brown & Rho can reproduce the new FOPI data in the 
covariant kaon dynamics quite well and the strength of the potential is roughly consistent 
with recent information obtained from p + A scattering [22], 


II. THE MODEL 


Due to its relativistic origin, the kaon mean field has a typical relativistic scalar-vector 
type structure. For the nucleons such a structure is well known from Quantum Hadron 
Dynamics [24], This decomposition of the mean field is most naturally expressed by an 
absorption of the scalar and vector parts into effective masses and momenta, respectively, 
leading to a formalism of quasi-free particles inside the nuclear medium [24], 

From the chiral Lagrangian the field equations for the iF -mesons are derived from the 
Euler-Lagrange equations [12] 


± + (< 


Drn 
p*2 P s 


0 K ± (x) = 0 


( 1 ) 


4 /* 2 ^“ ' V" /, 

Here the mean field approximation has already been applied. In Eq. (1) p s is the baryon 
scalar density, j f , the baryon four-vector current, f* the in-medium pion decay constant. 
Introducing the kaonic vector potential 
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Eq. (1) can be rewritten in the form [19] 


± iV^) 2 + mj? 0 K ±(^) = 0 . 


( 3 ) 


Thus, the vector held is introduced by minimal coupling into the Klein-Gordon equation. 
The effective mass of the kaon is then given by [3,19-21,25] 



Due to the bosonic character, the coupling of the scalar held to the mass term is no longer 
linear as for the baryons but quadratic and contains an additional contribution originating 
from the vector held. The effective quasi-particle mass defined by Eq. (4) is a Lorentz scalar 
and is equal for K + and K . In nuclear matter at rest the spatial components of the vector 
potential vanish, i.e. V = 0, and Eqs. (3) reduce to the expression already given in Ref. 
[12]. However, Eqs. (3) generally account for the correct Lorentz properties which are not 
obvious from the standard treatment of the kaon mean held [12-14,6]. 

The covariant equations of motion are obtained in the classical (testparticlc) limit from 
the relativistic transport equation for the kaons which can be derived from Eqs. (3). They 
are analogous to the corresponding relativistic equations for baryons and read [19] 


d jf = d J£l = + g , m 

dr mjk ’ dr rri^ 


( 5 ) 


Here = (t, q) are the coordinates in Minkowski space and = d^V u — is the 

held strength tensor for K + . For K~ where the vector held changes sign. The equation of 
motion are identical, however, F^ has to be replaced by —E 7 "'. The structure of Eqs. (5) 
may become more transparent considering only the spatial components 


dk* rnt drrit dV° k* ( d 

_ —____ _ -j- _ X I _ X y 

dt E* <9q <9q E* yclq 


( 6 ) 


where the upper (lower) signs refer to K + ( K~). The term proportional to the spatial 
component of the vector potential gives rise to a momentum dependence which can be 
attributed to a Lorentz force, i.e. the last term in Eq. (6). Such a velocity dependent 
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(v = k */E*) Lorentz force is a genuine feature of relativistic dynamics as soon as a vector 
field is involved. 

If the equations of motion are, however, derived from a static potential 
U(p, k) = u(p, k) — cu 0 (k) 

= ^k 2 + ml - j§TPs + V$±V 0 - ^/k 2 + ml (7) 

as given in Refs. [6,12-14,18], the Lorentz-force (LF) like contribution is missing. Non¬ 
covariant treatments are formulated in terms of canonical momenta k instead of kinetic 

momenta k* and then the equations of motion read 

dk rn^ dm ^ <917° <9V ?: 

dt E* <9q <9q * <9q 

with v = k */E* the kaon velocity. 

Following Brown and Rho [2], we use E K n = 450 MeV, f* 2 = 0.6 f 2 for the vector field 
and /* 2 = f 2 for the scalar part given by — Ekn/Z^Vs- Tins accounts for the fact that 
the enhancement of the scalar part using f* 2 is compensated by higher order corrections in 
the chiral expansion [2,4], This parameterization is denoted as Brown & Rho parameteri¬ 
zation (BRP), which has already been used in our previous investigations [19-21,25]. For 
comparison also a weaker potential with E K n=350 MeV and f* 2 = f 2 is applied. This pa¬ 
rameterization was originally used in ref. [12] and is called the Ko and Li parameterization 
(KLP) in the following. The baryon dynamics are treated within the framework of Quan¬ 
tum Molecular Dynamics (QMD). For the nuclear forces we use the standard momentum 
dependent Skyrme interactions corresponding to a soft/hard EOS (K=200/380 MeV). For 
the determination of the kaon mean field we adopt the corresponding covariant scalar-vector 
description of the non-linear cru model. Here we use the parameterizations of [13,26] which 
correspond to identical soft/hard nuclear EOSs. 

The potentials given in Eq. (7) at zero momentum by the BRP and KLP are shown in 
Fig. 1, where one also shows the kaon potential determined from the kaon-nucleon scattering 
length using the impulse approximation (IA) [27]. In IA the energy u(p, k), Eq. (7), is given 
as follows: 
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u(p,k) 


( 9 ) 



where rn^ is the nucleon mass and okn ~ -0.255 fm is the isospin-averaged kaon-nucleon 
scattering length in free space [28]. It is seen from Fig. 1 that the potentials predicted by 
BRP with the soft equation of state (solid curve) and by KLP with soft EOS (dashed line) 
are slightly stronger than that ones with the hard EOS (see the dot curve and long-dashed 
line, respectively). Up to saturation density the impulse approximation and the BR potential 
almost coincide (Uk(Po) — 30 MeV) but at supra-normal densities the BR potential rises 
much steeper than the IA. The KL potential (Uk(po) — 5 MeV, dashed line), on the other 
hand, is much weaker than those given by both, BRP and IA. 

The K + creation mechanism is treated as described in refs. [20,21,25] where one uses the 
improved cross section of refs. [29,30] for the baryon induced K + creation channels and of 
refs. [31] for the pion induced channels. The kaon production is treated perturbatively and 
does generally not affect the reaction dynamics [13,20,21,25,26]. The shift of the production 
thresholds of the kaons by the in-medium potentials are taken into account as described in 
[20,21,25]. In order to account for energy-momentum conservation it is useful to formulate 
the mass shell condition, Eq. (3), in terms of the canonical momenta 

0 = kg - mg = kl~ m 2 K - 2 m K U opt (10) 

with 

C V(,, k ) = -E s + T-M- + M. (ii) 

Here we introduced the total scalar kaon self-energy Eg = rn* K — mx- Since U opt is a Lorentz 
scalar it can also be absorbed into an effective mass 

rh K (p, k) = yjm& + 2m K U opt (p, k) (12) 

which sets the canonical momenta on the mass-shell 

0 = kg- mg = kl - t Ax • (13) 
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By definition rhy is a scalar but in contrast to who’s analog is the Dirac mass in the 
case of nucleons % absorbs the full optical potential and corresponds at zero momentum 
to the energy to. The threshold condition for K + production in baryon induced reactions 
reads then 


\fs > m B + my + rliK (14) 

with a/s the center-of-mass energy of the colliding baryons. The momenta of the outgoing 
particles are distributed according to the 3-body phase space 

rh B ,rh Y ,m K ) — d$ 2 (\/s, rh B , M)dM 2 $ 2 (M, rhy, fh K ) ■ (15) 


The two-body phase space in Eq.(15) has the form 

Trp*(y/s,mi,m 2 ) 


$2(Vs,mi,m 2 ) = 




(16) 


where 


p*(yfs,m 1 ,m 2 ) = 


y /(s — (mi + m 2 ) 2 )(s - (mi - m 2 ) 2 ) 


2-v/i 


(IT) 


is the momentum of the particles 1 and 2 in the c. m. frame. Eq. (15) corresponds to a 
distribution of the particle momenta according to an isotropic 3-body phase space. However, 
in [6] a parameterization of the form 


3 (v / s, 'm B) rh Y , m K ) = dW K (Vs, rh B , rhy, M K )dM^ 2 {\fs - M K , m Y , m B ) (18) 


has been suggested where the kaon momentum p is distributed according to 


dW K ~ 


P 


1 - 


P 


JP 


Pn 


(19) 


with p r n „ = p*(^/s, fh B + myprfiK) the maximal kaon momentum in the BB c.m. frame. 
Mk = y/p 2 + rh 2 K in Eq. (18). The parameterization of Eq. (18) has been motivated 
by analyzing corresponding pp pAK + data [6] and shifts the kaon spectrum to lower 
momenta compared to an ideal 3-body phase space. The optical potentials of the baryons 
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which enter via rfiB, my are taken from the soft/hard EOS versions of the au model [13,26]. 
The hyperon fields are thereby scaled according to SU(3) symmetry U° pt = | U^ t . Since the 
m’s depend on the final state momenta the determination of d<f> 3 is a self-consistency problem 
which is solved by iteration. The same procedure is applied to the 2-body-phase-space in 
pion induced reactions. The rescattering of the K + mesons with baryons and the Coulomb 
interaction are taken into account. The electromagnetic interaction is treated analogously 
to the strong interaction, i.e. by adding F^ u = d^A 1 ' — d u A 11 given by the electromagnetic 
vector potential to the Eq. (5). 


III. RESULTS 

In order to study first the influence of the 3-body phase space in BB collisions we 
consider in Fig. 2 C + C collisions at 2.0 AGeV. The KaoS Collaboration has measured 
inclusive K + spectra at various c.m. angles with high precision. Corresponding calculations 
are performed with the soft EOS including the K + in-medium potential BRP and final 
state interactions by rescattering processes. In this scenario the excitation function of the 
kaon multiplicities measured by KaoS [10] is well reproduced over a large energy range from 
0.8 to 2 AGeV [20]. Without in-medium potentials the total multiplicities are significantly 
overpredicted. Looking to the spectra in more detail we find that an isotropic 3-body 
phase space in the BB —> BYK + channel shifts the spectra to too high momenta. The 
empirical parameterization of the 3-body phase space suggested in ref. [6], Eq. (17), improves 
the situation but is still not fully sufficient to account for the angular asymmetry which 
is seen in the KaoS data [34], A relatively good agreement can be achieved introducing 
an empirical angular dependence da oc (1 + a cos 2 6 c , m )dcos 6 c . m . in the elementary cross 
sections. Following the data analysis of [34] an asymmetry parameter a = 1.2 leads to slightly 
forward/backward peaked elementary BB —> BYK + cross sections and the corresponding 
spectra shown in Fig. 2 are then well reproduced by the transport calculations. 

In order to investigate the influence of covariant dynamics on the K + in-plane flow we 



consider the 1.93 AGeV 58 Ni + 58 Ni collisions at impact parameter b < 4 fin, corresponding 
to the FOPI centrality cut. The K + rapidity distributions dN/dY 0 (Y° = Yi ab /Y C M -1) for 
this reaction is shown in Fig.3. Again the BRP in-medium kaon potential has been used. 
Consistent with our previous results [20] and those from other groups [16], this figure shows 
that the existence of a repulsive kaon potential is required to match the experimental data 
[9,11], in particular around mid-rapidity. The influence of the nuclear EOS is thereby on the 
20% level. The best agreement with the data is obtained with the soft EOS and including 
a kaon potential. 

The next figure gives the corresponding proton flow. A p t cut of p t /m > 0.5 accounts 
for the experimental acceptance [7,23]. At this high energy the dependence of the proton 
flow on the nuclear EOS is weak and both options, i.e. a soft and a hard EOS reproduce 
the FOPI data [7] reasonably well. It should be noted that the analysis has been performed 
for all protons. At large rapidities the experimental flow is rather well reproduced whereas 
deviations occure between — 1 < Y° < 0 where the consideration of cluster effects might 
help to improve on the data fitting. It is, on the other hand, known from Ref. [15] that 
proton and kaon flow are only very loosely connected. The kaon producing sources carry 
this large in-plane flow. The kaons themselves carry a much smaller flow fraction since, at 
given rapidity, they are produced from two baryons originating from very different rapidity 
regions. Therefore baryon sources with positive and negative p x /m values add up to an 
almost vanishing net flow (see the results given by open down triangles in Fig.5). A’s 
which are produced in association with Jl + , s have very similar flow pattern as protons, i.e. 
they show almost the same p x /m scaling [23]. In [32] the A flow was investigated within 
the present model and the data were best reproduced including the hyperon mean field 
according to SU(3) scaling Uj t = | U^ pt . As also observed in [33] the primordial A flow is 
moderate but strongly enhanced by the AN final state interactions, i.e. rescattering and the 
mean field. 

Figure 6 shows now the K + transverse flow as a function of the scaled rapidity in 1.93 
AGeV 58 Ni + 58 Ni reactions. In this figure the full squares represent the ’95 data set from 
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FOPI (old data) [7]. The full circles show the ’99 data from FOPI [23] with improved 
statistics, open circles indicate their reflection at mid-rapidity. The theoretical results are 
given by the BRP with a soft EOS. Around mid-rapidity the two calculations with [/#& LF 
(static potential + Lorentz force, open down triangles) and without Uk (open up triangles) 
almost coincide. The result without Uk, he. only including the kaon rescattering effect, 
predicts a slightly positive flow. The result with Uk& LF leads to a very small anti-flow. 
Around mid-rapidity both calculations agree with the data within error bars. However, 
at spectator rapidities the two results with Uk& LF and without Uk differ substantially 
from each other. With respect to the old data set both calculations, i.e. with Uk& LF 
and without Uk, agree with experiment within error bars since both reproduce the nearly 
vanishing side flow signal of Jl + s around mid-rapidity. This means that the old data can be 
reproduced without need for an in-medium kaon potential [15]. However, the new data with 
much small error bars are only in agreement with those calculations which treat the kaon 
potential in the covariant kaon dynamics. This indicates that it is necessary to include the 
in-medium kaon potential in order to describe the new data for the K + transverse flow. 

From Fig. 6 it is seen that the strongly repulsive static potential tends to push the 
kaons dramatically away from the spectator matter, leading to a strong anti-flow around 
midrapidity (stars). The effect of the LF contribution in the covariant kaon dynamics pulls 
the kaons back to the spectator matter, resulting in a finally reasonable pattern of the K + 
transverse flow, which is in good agreement with the FOPI data. This feature of the LF 
contribution can also be seen in Fig. 7 where the calculations are performed by the BRP 
with a hard EOS. This illustrates that the LF like contribution, originating from spatial 
components of the vector field, provides an important contribution to the in-medium kaon 
dynamics in heavy ion collisions. Kaons are produced in the early phase of the reaction 
where the relative velocity of projectile and target matter is large. Thus the kaons feel a 
non-vanishing baryon current in the spectator region, in particular in non-central collisions. 
This contribution dramatically counterbalances the influence of the repulsive potential on 
the Ii + transverse flow. 
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From the comparison between the results in Figs. 6 and 7, on the other hand, it is found 
that the K + transverse flow calculated with a hard EOS gives rise to a slightly positive flow 
around mid-rapidity and a worse fit to the new data in the target region. This is due to 
the fact that the kaon in-medium potential calculated with a hard EOS is weaker than that 
one given with a soft EOS (see Fig.l). However, the K + transverse flow given with a hard 
EOS can also roughly describe the new data set within error bars. Thus it appears to be 
impossible to extract information on the nuclear EOS from the analysis of the K + side flow. 

The K + transverse flow given by the KLP with a soft EOS is shown in Fig. 8 where 
the results by the BRP with Uk& LF for the soft EOS (open down triangles) are shown as 
well. In this figure the experimental data and their reflections with respect to mid-rapidity 
are represented by the same symbols as in Fig. 6. The results (stars) are calculated by 
the KLP with Uk but without (w/o) LF. The standard treatment with a static potential, 
Eq. (7), shows an anti-flow around mid-rapidity which is larger than that one given by the 
BRP with Uk& LF and is in better agreement with the old data. The situation changes, 
however, dramatically when the full Lorentz structure of the mean held is taken into account 
(diamonds). The influence of the repulsive potential of K + s on the in-plane flow is almost 
completely counterbalanced by the velocity dependent part of the interaction, i.e. the LF 
contribution. Hence, no net effect of the potential is any more visuable. This is clearly 
seen from the comparison between this result (diamonds) and the ones (open up triangles) 
obtained w/o Uk- This seems to indicate that the static potential predicted by the KLP is 
too weak to overcome the cancellation effects of the LF on the flow. Using the BRP, the 
kaon in-medium potential has a stronger repulsive character which gives rise to a reasonably 
anti-correlated flow signal and leads to a good agreement with both, the old and new data 
set from FOPI. 

In this context one has to keep in mind that the net source flow for the kaon production 
is experimentally not accessible. Here one has to rely on the predictions from transport 
models. The analysis of the small kaon flow signal which arises to large amount from the 
cancellation of large source flow components carries inherent uncertainties due to experimen- 
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tal and theoretical error bars of the source flow. Investigating the influence of the in-medium 
kaon potential we were confronted with two effects: If the kaon dynamics is described by 
the standard treatment with a static potential, i. e. neglecting the effect of a Lorentz force 
from spatial component of the vector held, the calculated K + in-plane how is not able to 
reproduce the new FOPI data with reasonable precision. The Lorentz force effect is large 
and it changes the K + how patter qualitatively. We consider this qualitative result therefore 
as stable. Quantitative statement on the strength of the kaon potential have to be made 
more carefully. Since differences are most pronounced at large spectator rapidities where 
the experimental proton how is within error bars well reproduced by the present transport 
calculations, the K + how can help to distinguish between different models. However, clue 
to the smallness of the signal and present theoretical and experimental error bars, def ini te 
quantitative statements require further investigations. 

IV. CONCLUSIONS 

In summary, the influence of the chiral mean held on the transverse how of kaons in 
heavy ion collisions at SIS energy has been investigated within covariant dynamics. The 
kaons inside the nuclear medium are described as quasi-particles carrying effective masses 
and momenta. This accounts for the correct mass-shell properties of the particles inside 
the nuclear medium. A consequence of the covariant kaon dynamics is the appearance of 
a momentum dependent force proportional to the spatial components of the vector held 
which resembles the Lorentz force in electrodynamics. Although such Lorentz forces vanish 
in equilibrated nuclear matter they provide an essential contribution to the dynamics in the 
case of energetic heavy ion collisions. The inhuence of in-medium effects on the in-plane 
how is counterbalanced to large extent by this additional contribution. 

In the present studies two types of in-medium potentials with different parameterizations 
have been applied. Their inhuence on the K + in-plane how has been discussed. Our theo¬ 
retical results show that in the covariant dynamics the new FOPI data can be reasonably 
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described using a parameterization proposed by Brown & Rho which partially accounts for 
higher order corrections in the chiral expansion. The reproduction of the more recent FOPI 
data, in particular at spectator rapidities, requires a relatively strong repulsive K + poten¬ 
tial, which is in good agreement with the ones determined from the kaon-nucleon scattering 
length using the impulse approximation [27] and also with recent information obtained from 
p + A reactions [22]. This finding is consistent with the description of the K + multiplicities 
in A + A reactions. Most transport simulations reproduce corresponding data [8] only when 
in-medium potentials are included [12,14,16,20]. In summary, information extracted from 
the analysis of the K + transverse flow is consistent with the knowledge from other sources. 
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FIGURES 



FIG. 1. Density dependence of the in-medium kaon potential at zero momentum. The KLP 
and BRP stand for the Ko&Li parameterization and Brown&Rho parameterization, respectively. 


17 







_ I _I_I_I_I_I_I_I_ii_I_ 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

P lab [GeV/c] 


FIG. 2. Inclusive K + spectra at 9 i a b = 32°, 40° in 12 C + 12 C reactions at 2.0 AGeV. The 
calculations are performed using different descriptions of the BB —► BYK + final state 3-body 
phase space. The dotted curves refer to an isotropic 3-body phase space while the solid curves 
are obtained using the parameterization of Eq. (19) with an additional empirical c.m. angular 
anisotropy (see text). The corresponding KaoS data are taken from [34], 
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FIG. 3. K + rapidity distributions in 1.93AGeV 58 Ni + 58 Ni reactions at impact parameter 
b< 4 fm. The calculations are performed with and without the BRP K + in-nredium potential and 
using a soft/hard EOS. Diamonds represent data from KaoS [9], triangles those from FOPI [11]. 
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Y (0) 

FIG. 4. Transverse proton flow as a function of rapidity Y° in 1.93AGeV 58 Ni + 58 Ni reactions 
at impact parameter b< 4 fm. The calculations with a soft (y) and hard (A) nuclear EOS are 
compared to the FOPI data [7]. 
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FIG. 5. Transverse flow of the kaon producing sources as a function of rapidity Y° in 1.93AGeV 
58 Ni + 58 Ni reactions at impact parameter b< 4 fm. The calculations are performed with the BRP 
K + in-medium potential and using a hard nuclear EOS. The FOPI data are from [7]. 
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FIG. 6. The K + transverse flow as a function of rapidity Y° in 1.93AGeV 58 Ni + 58 Ni reactions 
at impact parameter b< 4 fm. The calculations are performed using the BRP K + in-medium 
potential with a soft EOS. The full squares represent the FOPI data from [7], full circles are 
more recent FOPI data [23], their reflections with respect to mid-rapidity are shown by the open 
circles. The open down triangles denote the calculated results with Uk and the Lorentz-force (LF) 
contribution. The open up triangles indicate the results without Uk- The stars stand for the 
results with Uk but without LF. 
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FIG. 7. Same as Fig. 6, but with a hard EOS. 
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Uk& LF are also shown. 
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